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1 Introduction 


Last decade was very fruitful in spectroscopy of the rapidly growing number of particles. 
Many new particles are discovered (see [1-9] and the references therein) which part of 
them are successfully described in framework of the quark model, and interpretation of 
the remaining ones demand going beyond the quark picture. In particular, the masses and 
decay widths of the excited mesons such as ©2(2400), ©5^(2573), and 13^2(5840) 

with the quantum number = 2^ are measured in experiments [8-10]. The properties of 
©-wave and radial excited heavy light meson system will be examined in detail in future 
planned experiments at LHBc and KEK-B. 

The strong decays of ©-wave mesons, such as ©2(2460) —©*'''7r“, ©+7r“, ©^(2460) 
© 07 r+ [6,12-14], ©1^(2573) ^ ©°iL+ [6], 132(5747) ^ i3*+7r-, iS+Tf, [9,10], and i3+(5840) ^ 

K~ [9,10] have already been observed in experiments. The decay constants of the tensor 
mesons ©^(2460) and ©1^(2573) have been studied within the three-point QCD sum rules 
method in [15]. In framework of the same approach the strong constants of ©2(2460) —)■ ©tt, 
©^^(2573) —)■ ©77, 132 ( 5747 ) —)■ Bn, and 13^2(5840) —)■ BK transitions have been calculated 
in [16]. Recently, the strong decay constants of the ©^(2460) —)■ ©*7r, ©^2(2573) —)■ ©*77, 
132 ( 5747 ) —)■ 13* TT, and 135^(5840) —)■ B*K transitions have been studied in [17] in framework 
of the three-point QCD sum rules and the local QCD methods. 

The present work is devoted to the calculation of the constants of the tensor mesons 
© 2 ( 2460 ), ©5^(2573), 133 ( 5747 ), and 13^2(5840) in framework of the QCD sum rules method. 

The paper is arranges as follows. In section 2, the light cone QCD sum rules derived 
for the coupling constants of the ©2 —)■ ©tt, ©^^ —t ©77, B^ —)■ Bn, and Bg^ —)■ BK 
transitions. Section 2 is devoted to the numerical analysis of the aforementioned coupling 
constants, and present the values of the corresponding two-body strong decays. In this 
section, comparison of our results with the ones existing literature, and our concluding 
remarks are also presented. 


2 Light cone QCD sum rules for the strong coupling 
constants of heavy tensor mesons 

In this section we derive the light cone QCD sum rules for the strong coupling constants 
of heavy tensor mesons with heavy pseudoscalar (vector) mesons with participation of the 
light pseudoscalar mesons. For this aim we consider the following correlation function. 
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where 


= o Vu Q{x) + q{x)')y Q{x) 


is the interpolating current of the heavy tensor meson, J 5 = qi'y^Q, and © = q'jrQ are 
the interpolating currents of the heavy pseudoscalar and vector mesons, respectively; © 
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denotes the light pseudoscalar meson, q and Q in the expressions of the interpolating 
currents describes the light and heavy quarks. The covariant derivative is dehned in 
the following way, 
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where A“ are the Gell-Mann matrices, and is the external gluon held, and in our 

calculations we shall use the Fock-Schwinger gauge, i.e., x^A^^{x) = 0. The advantage of 
this gauge is that the gluon held can be expressed in terms of gluon held strength tensor 
dehne as. 


Al= [ dtGuf,(tx)x'' 

Jo 

= + ^x^x'^V^G,^ + • • ■ . (4) 

In order to construct the sum rules for the strong coupling constants of the heavy tensor 
mesons with heavy pseudoscalar (vector) mesons with the participation light pseudoscalar 
mesons, we shall calculate the correlation function (1) in two diherent ways: a) In terms of 
hadronic parameters; b) in terms of quarks and gluons. Equating then these two represen¬ 
tations, we derive the sum rule for the aforementioned hadronic coupling constant. 

We hrst consider calculation of the correlator function from the hadronic side. Inserting 
a complete set of hadrons carrying the same quantum numbers as the interpolating currents 
and J 5 (t-) into the correlation function, and isolating the ground state contributions of 
the heavy tensor and heavy pseudoscalar (vector) mesons we get. 
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where Tq, Pq,Vq,V mean heavy tensor, heavy pseudoscalar, heavy vector mesons, and light 
pseudoscalar mesons, respectively, and dots mean higher states and continuum contribu¬ 
tions. 

The matrix elements in Eq. (5) are dehned as. 
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where is the polarization vector of the heavy tensor meson, fx, fq, fv are the decay 
constants of the tensor, heavy pseudoscalar, and heavy vector mesons; gi and g 2 are the 
hadronic coupling constants, and s is the polarization vector of the heavy vector mesons. 
Performing summation over the spins of the heavy tensor and heavy vector mesons in Eqs. 
(5-10), we get for the correlation functions describing the heavy tensor mesons with heavy 
pseudoscalar, heavy vector mesons and light pseudoscalar mesons, the following results, 
respectively. 
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where 
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We now turn our attention to the calculation of the correlation function from the QCD 
side. i.e. in terms of the quark-gluon degrees of freedom. After contracting the heavy 
quark helds the correlation function can be written as. 
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where Sq{x) is the heavy quark propagator whose expression in the coordinate representa¬ 
tion is given as. 
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where Kn is the modihed Bessel function. In further analysis we shall use the following 
relations for the derivatives of the Bessel functions, 

d 1 r 

dx 2 L J 

^{x^Kn) = -x-Kn-i{x) , 

^{x-^K^) = -x-^K,,+,{x) . (15) 

dx 

It should be noted here that the expansion of the quark operator up to twist four terms 
is calculated in [18], which gets contributions from qGq, qGGq, and qqqq nonlocal operators. 
In the present work we consider operators with one glnon held. Neglecting the contribntions 
coming from the qGGq and qqqq operators can be justihed on the basis of an expansion in 
conformal spin. 

Having the expression of the heavy qnark operator at hand, we now calculate the cor¬ 
relation fnnction from the QCD side. The expression of the correlation fnnction in deep 
Encledian domain —)■ —oo, p'^ —)■ —oo, can be obtained by nsing the operator prodnct 
expansion. It follows from Eq. (13) that in order to calculate the correlation fnnction 
from the QCD side the matrix elements of the nonlocal operators between vacnnm and 
one-particle light psendoscalar meson states. The matrix elements {V{q) |g(a;)rg'( 0 )| 0) are 
parametrized in terms of the distribntion amplitndes (DAs) [19, 20] as given below. 
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and qi and q 2 are the quarks in the meson V, Va = daqdaqdagSil — aq — aq — Ug), and 
and the DA’s A(m), M{u), ^riu), 0p(m), 0o-(m), Ticii), x4,_L(Q;j), ^ 11 ( 0 ;*), V±iai) and V||(Q;j) 
are functions of dehnite twist and their expressions are given in the next section. 

Using Eq. (16) and corresponding coefficients of the Lorentz structures the theoretical 
part of the correlation function can be calculated straightforwardly. Equating the coeffi¬ 
cients of the respective Lorentz structures calculated from the hadronic and QCD sides of 
the correlation function, and performing Borel transformation for the variables —p^ and 
—p'2 = _|_ g )2 order to suppress the continuum and higher states contribution, we 

obtain the sum rules for the strong coupling constants of the heavy tensor with pseudoscalar 
heavy and light mesons given as, 
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where~means Borel transformed invariant function, Pq and Vq refer to the heavy pseu¬ 
doscalar and vector mesons, respectively. The expressions of Ai and Bi in the above¬ 
expression are as follows, 
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1, for the e^^^^'^PaquQr] structure , 
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The expressions of the and are given as: 
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5) Coefficient of the PaQu^ri and e’'°‘'^'^paqf_iqri structures 
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The function ji(/(«)) appearing in the invariant functions above is defined as: 
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It should be noted here that, in the above expressions the light quark masses rriu, nid 
and rus are all set to zero, while in the numerical calculations the mass mg of the strange 
quark is taken into account. It should also be mentioned that the derivation of the double 
spectral density for the higher twist contributions is not trivial, and it is calculated in [21], 
In our calculations the higher twist contributions appear in the similar form as they do 
in [21]. Therefore, for more detail about this issue, the interested reader can consult this 
work. 


3 Numerical analysis 

This section is devoted to the numerical analysis of the sum rules for the heavy tensor 
mesons with the heavy pseudoscalar (vector) and light pseudoscalar mesons. 

The main input parameters of the light cone QCD sum rules are the distribution am¬ 
plitudes (DAs), whose expressions are given below [19,20], 
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The valnes of the parameters of, of, 7 / 3 , 7 / 4 , w^, and 104 entering Eq. (20) 
Table (1) for the psendoscalar n, K and 7 / mesons. 
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a[ 

0 

0.050 

02 (set-1) 

0.11 

0.15 

02 (set-2) 

0.25 

0.27 

V3 

0.015 

0.015 

h4 

10 

0.6 

W3 

-3 

-3 

W4 

0.2 

0.2 


Table 1: Parameters of the wave function calculated at the renormalization scale /i = 1 GeV 

In addition to the above-mentioned parameters, the decay constants of the heavy tensor, 
heavy pseudoscalar (vector) and light pseudoscalar mesons, and masses of the quarks. The 
decay constants of the heavy tensor mesons are calculated in [16,17,22] to have the values 
fv* = (0.018 ± 0.007), = (0.023 ± 0.011), h* = 0.011, = 0.013, = 0.24 Gef/, 

/b* = 0.16 GeV" [21]. ' _ 

In the present work we use the MS values of the quark masses predicted by the 
particle data group whose values are [8]: fficirnc) = (1.275 ± 0.025) GeV", fnh{mb) = 
(4.18±0.03) GeV", and 771^(2 GeV") = (0.095±0.005) GeV". The masses of the heavy mesons 
we use in the present work are calculated in framework of the the QCD sum rules, having the 
values, mw = (2.460 ± 0.009) GeV", mm = (5.73 ± 0.06) GeV", = (5.84 ± 0.06) GeV", 
which are all very close to their experimental values. 

There are two extra parameters entering to the sum rules, namely, the continuum 
threshold sq and the Borel mass parameters and M|. In the present analysis we set 
Mf = M| = 2M^, and this choice corresponds to uq = 1/2. The continuum threshold sq is 
determined from an analysis of the two-point correlation function which leads to the follow¬ 
ing results: So,x>* = (8.5 ±0.5) GeV"^, = (9.5 ±0.5) GeV"^, = (39.0 ±1.0) GeV"^, 
= (41.0 ± 1.0) GeV"2 [16,17]. 

Of course, we need to hnd such a region of where the strong coupling constants 
9TqPqV 9TqVqV are insensitive to the variation in M^. The upper bound of M^is deter¬ 
mined from the condition that the higher states and continuum contributions constitute, 
say, 30% of the total result. The lower bound of the Borel mass parameter is deter¬ 
mined from the condition that the contribution of the highest term with the power 1/M^ 
is less than 25% of the contribution coming from from the highest power of M^. These two 
conditions lead to the following working regions of 2 GeV"^ < < 4 GeV"^ (for 

and T>sJ, 4 GeV^ <M^ <7 GeV^ (for and SJJ. 

Having decided the working regions M^, we now calculate the strong coupling constants 
for the TqPqV and TqVqV vertices. Performing similar analysis we obtain the values of 
the coupling constants for the TqPq{Vq)V vertices which are presented in Tables 2 and 3. 
In these calculations we use two different sets of parameters that appear in the expressions 
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of the DAs. 
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18.3 
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12.2 
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78 ±20 

74 ±5 

10 ± 1 

10 ± 1 

0.10 

0.13 

19 ±6 

9±3 


- 

- 

39.3 

17.1 

90 ±25 

88 ±20 

70± 15 

73 ±20 

0.30 

0.15 

40± 12 

36± 12 

Bs2 BK° 

- 

- 

26.3 

12.9 

400 ± 100 

360 ± 50 

58± 12 

55± 12 

0.25 

0.30 

24 ±8 

22 ±8 


Table 2: The values of the strong coupling constants for the TqPqV vertices for different 
Lorentz structures 
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73 
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4.00 

0.73 

4.4 ± 1.0 
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0.42 ±0.10 

no stability 

^ V*-K+ 

2.98 

0.79 

4.4 ± 1.0 

no stability 

-0.37 ±0.13 

no stability 


3.87 

- 

5.9 ±2.0 

no stability 

0.22 ±0.05 

no stability 

Bl^ B*+K- 

2.89 

- 

3.0 ± 1.0 

no stability 

0.28 ±0.06 

no stability 


Table 3: Same as in Table 2, but for the TqVqV vertices 

In Table 3 “no stability” means that there is no region of the Borel parameter where 
the results for the strong coupling constants are insensitive to its variation. 

From these Tables we deduce the following conclusions: 

• The value of the strong coupling constant depends very strongly on the choice of the 
corresponding Lorentz structure, especially for the Tq —)■ PqV transition. 

• From our analysis we observe that the value of the strong coupling constant for the 

Tq —)■ PqV transition ranges in a rather wide region from 0.3(0.15) to 290(400) for 
the 82 —^ —)■ In the case of T> 2 (S 2 ) (V^K^) transitions the 

valnes of the corresponding coupling constants vary in the range 0.15 to 78.0. These 
resnlts also show that our analysis that the valnes of the strong conpling constants 
are very sensitive to the valnes of the parameters appearing in the DAs. 

It follows from Table 1 that the most reliable valne for the TqPqV vertex follows from 
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the q^pu + QuPfi structure, from which we get, 


{ 17 ±5, ^ V+n- , 

19 ± 6, ©+ ^ V^K+ , 
40 ± 12, ^ i3-7r+ , 

24 ± 8, B^K^ . 


( 21 ) 


The strong coupling constants for the Tq —)■ VqV transitions is obtained for the EpaarP^p'^q'^ 
structure, and for the Set 1 values of the parameters of the wave functions are used. As 
the result, for the coupling constants of the TqVqV vertex we get. 


{ 4.4 ±1.0, Vl^V*+Ti- , 
4.4 ± 1.0, T>+ ^ V*^K+ , 
5.9 ± 2.0, Bl ^ B*-ti+ , 
3.0 ± 1.0, Bl^ -> B*+K- . 


( 22 ) 


We now compare our results on the strong coupling constants with those predicted by 
the 3-point QCD sum rules method. For the structure our prediction for the coupling 
constant gv^v+ir- is approximately two times larger compared to the one predicted in [15] 
and [16]. For the B 2 —)■ B~n~^ and Bs 2 —t transition coupling constants our results 

are two and more than ten times larger than the ones predicted in [16]. 

At the end of this section we present the results for the decay widths of all above- 
considered transitions, whose expressions are gives as. 


. gl bf 
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where 
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Note that the |p|^ dependence is an indication of the fact that the decay takes place at the 
D-wave level. The results for the decay widths can be snmmarized as follows. 


• In the case for the transitions of the heavy tensor mesons to heavy psendoscalar and 
light psendoscalar mesons (for the p^q^ + P^q^ strnctnre). 


F (T'°(2460) ^ P+TT-) = 8.6 x 10"=^ GeV , 

F (T>+(2573) ^ = 4.4 x 10"=^ GeV , 

F (B^(5747) ^ i3-7r+) = 3.7 x lO'^ GeV , 

F (13^^(5840) ^ B^K^) = 8.6 x lO"® GeV . (25) 
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• In the case for the transitions of the heavy tensor mesons to heavy vector and light 
pseudoscalar mesons (for the structure), 


r (r'°(2460) ^ V*^Ti-) = 5.2 X 10"^ GeV , 

r (r>+(2573) ^ = 2.3 x 10"^ GeV , 

r (^ 2 ( 5747 ) ^ i3*-7r+) = 4.0 X 10"^ GeV" , 

r (13^^(5840) ^ = 6.9 x 10“^ GeV . 


It follows from the experimental data that [8], 

r(p2(2460) ^ V+n- 


r(r>2(2460) ^ V*+Ti-) 
r(i32(5747) ^ g+TT-) 

TT-) 


= 1.55 , 
= 0.91 , 


r(i32(5747) ^ B 

as well as from the BaBar Collaboration data [23] 
r(p2(2460) ^ P+TT-) 

r(r>2(2460) ^ V+Ti-) + r(r>2(2460) ^ V*+n-) 

When we calculate the same ratios from Eqs. (25) and (26), we obtain that, 

r(r>2(2460) ^ V+n- 


= 0.62 ±0.03 ±0.02 . 


and 


r(r>2(2460) -> V*+Ti- 
r(i32(5747) ^ i3-7r+ 
fOB^(5747)“H^S^ 

r(p2(2460) -> P+TT") 


= 1.64 , 
= 0.93 , 


r(r>2(2460) ^ V+n-) ± r(r>2(2460) ^ V*+ti-) 


= 0.67 . 


We see that our prediction on the ratio of the decay widths are in good agreement with 
the experimental results, as well as they are quite close to the results of the works [16] and 
[17], 

Few words about the perturbative 0{as) corrections are in order. These corrections 
increase the correlation function of the coupling constant of the B* —)■ Bn transition about 
50% in the light cone QCD sum rules [24]. If we assume that this increase in the correlation 
function is correct, then the we expect that the coupling constant increases at the same 
order, i.e., gi —)■ l.5gi. This increase in the coupling constant doubles the values of the 
decay widths as well. 

In conclusion, we calculate the strong coupling constant of the heavy tensor meson to 
the heavy pseudoscalar (vector) and light pseudoscalar mesons. It is seen that the values 
of the hadronic decay constants are very strongly dependent on the choice of the Lorentz 
structure. Furthermore, using the value of the coupling constants calculated in this work 
we also estimate the corresponding decay widths. A comparison of our predictions on these 
hadronic coupling constants with the results of the 3-point sum rules is presented. 
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